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Here we present first results on a non-radiolytic, thermally-driven reaction sequence in solid H 2 0 + S0 2 + H 2 0 2 mix- 
tures at 50-130 K, which produces sulfate (SOj - ), and has an activation energy of 53 kj/mole. We suspect that these 
results may explain some of the observations related to the presence and distribution of H 2 0 2 across Europa’s surface 
as well as the lack of H 2 0 2 on Ganymede and Callisto. 
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1. Introduction 

The radiation-driven chemistry of surface and near-surface ices of Europa has 
been studied by several research groups in recent years. In contrast, thermally-dri- 
ven reactions in outer Solar System ices have remained relatively unexplored. In 
this note we report new laboratory experiments demonstrating that a straightfor- 
ward sequence of hydrolysis, protonation, and redox chemistry can convert S0 2 into 
SO 4 - (sulfate) at the temperatures of Europa’s surface and in the presence of H 2 0 2 
(hydrogen peroxide), a well-known radiation product of H 2 0-ice. These reactions 
suggest that sub-surface SO 4 - formation will occur in the absence of the direct 
reach of jovian magnetospheric radiation. 

Laboratory studies of H 2 0-ice, one of the more abundant substances on cold 
surfaces in outer space, have shown that radiation can cause both amorphization 
of crystalline H 2 0-ice (Hudson and Moore, 1995; Strazzulla et al., 1992) and the 
production of new molecules, such as H 2 0 2 , H 2 , and 0 2 . Of these three species, 
H 2 0 2 (hydrogen peroxide) formation has been observed in situ by laboratory 
infrared (IR) spectroscopy in H 2 0-ice irradiated at temperatures relevant to 
extraterrestrial surfaces (Loeffler et al., 2006a; Moore and Hudson, 2000; Zheng 
et al., 2006). Although all of these prior studies showed that the H 2 0 2 abundance 
is relatively low in irradiated H 2 0-ice, the H 2 0 2 absorption at 3.5 pm is well- 
separated from the stronger OH stretching vibrations of water, making it an 
attractive candidate for observation on planetary surfaces if both H 2 0-ice and 
radiation are present. 

Given the many extraterrestrial objects with surficial H 2 0-ice exposed to radi- 
ation, one might expect to find H 2 0 2 on a great number of icy bodies. However, 
H 2 0 2 has been detected unequivocally only on Europa (Carlson et al., 1999), with 
IR spectra of Ganymede and Callisto failing to show the 3.5-pm peroxide band. 
Interestingly, ultraviolet spectra of Ganymede do show 0 3 (Noll et al., 1996), a reli- 
able indicator of radiation processing of H 2 0-ice containing 0 2 . Differences exist 
among surface environments of Europa, Ganymede, and Callisto, such as concerning 
radiation flux, average surface temperature, and composition, but so far none of 
these factors readily explains why Europa would be the only place where H 2 0 2 is 
present. 
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In recent years we have used IR spectroscopy to investigate the solid-phase 
radiation chemistry of H 2 0 + S0 2 ices (Moore et al., 2007), solid sulfuric acid, 
and sulfuric acid hydrates (Loeffler and Hudson, 2012; Loeffler et al., 2011), all 
at temperatures relevant to the icy jovian satellites. With H 2 0 + S0 2 ices, it 
was observed that radiation readily produces SOj - (sulfate ion), presumably by 
the oxidizing action of radiolytically-generated H 2 0 2 . This assumption has been 
tested in the present paper by studying IR spectra of unirradiated H 2 0 + S0 2 + H 2 _ 
0 2 ices as they are warmed. We have observed that even at relatively low tem- 
peratures, H 2 0 2 and S0 2 react thermally in the presence of H 2 0-ice. The main 
products of the reaction have been identified and the reaction’s rate has been 
measured at 110-122 K, temperatures found on the jovian icy satellites men- 
tioned above. 

2. Experimental methods 

Experiments were performed with a cryostat (T m in ~ 10 K) operating in a stain- 
less steel high-vacuum chamber (P ~ 1 x 10 -7 Torr). Ice films were prepared by co- 
deposition of H 2 0, H 2 0 2 , and S0 2 onto a pre-cooled (50-100 K) gold-coated alumi- 
num mirror using three separate pre-calibrated gas lines. Pure H 2 0 2 was prepared 
in a glass manifold, using the technique described previously (Loeffler and Baragi- 
ola, 2011). During deposition, the sample’s thickness was monitored with interfer- 
ometry using a diode laser (670 nm). In all experiments the thickness of the 
resulting H 2 0 + S0 2 + H 2 0 2 (80;14:6) ice was 1.3 ±0.1 pm, assuming indices of 
refraction at 670 nm are similar to those measured in the visible region, 1.31 for 
H 2 0 (Merwin, 1930), 1.4 for H 2 0 2 (Giguere, 1943) and 1.36 for S0 2 (Musso et al., 
2000), uniform mixing, and an ice density of 0.96 g/cm 3 (0.82 g/cm 3 for H 2 0 (West- 
ley et al., 1998), 1.6 g/cm 3 for H 2 0 2 (Loeffler et al., 2006b) and 1.49 g/cm 3 for S0 2 
(Weast, 1983). 

Each sample’s IR spectrum was recorded before, during, and after warming at 
1 K/min to the annealing temperature (110-130 K). Spectra of ices were measured 
from 7000 to 400 cm -1 with a Bruker Vector Fourier Transform infrared spectrom- 
eter at 2-cm _1 resolution and with 200-scan accumulations. To obtain a spectrum, 
the reflectance from the ice-coated substrate was ratioed against the reflectance of 
the bare substrate, taken before ice formation, and then converted to absorbance 
units. 
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To quantify the kinetics of the thermal reactions occurring in our samples, we 
first fit the baseline of the 2850 cm -1 overtone band of H 2 0 2 with a non-linear curve 
and then integrated this same H 2 0 2 IR feature. We chose to examine H 2 0 2 and not 
S0 2 because the sublimation of pure H 2 0 2 is negligible below 130 K, whereas the 
sublimation of pure S0 2 can become significant above 100 K, the temperature re- 
gion of interest. We note that when H 2 0 2 is dispersed in water the strength of 
the 2850 cm -1 band depends slightly on temperature (Loeffler et al., 2006a). This 
dependence was confirmed in H 2 0-H 2 0 2 calibration experiments and we made 
use of it in our calculations; the strength of the H 2 0 2 band in question drops linearly 
by about 20% between 50 and 130 K. 

3. Results 

Fig. 1 shows the IR spectrum of an H 2 0 + S0 2 + H 2 0 2 (80:14:6) ice after deposi- 
tion at 50 K, during warming to 1 1 8 K, and while held at 1 1 8 K. At 50 K, each of the 
three molecules of the sample is easily identified in this spectral region: H 2 0 2 
(2845 cm -1 ), H 2 0 (1655 cm -1 ), and S0 2 (1329 and 1151 cm -1 ). On warming to 
118 K, new absorptions, belonging to HSOJ (1035 cm -1 ), SOj - (1092 cm -1 ), and 
S 2 Og - (954 cm -1 ) appeared, as did a broadening on the high-wavenumber side of 
the H 2 0 band at 1655 cm -1 , indicating H 3 0 + formation. As the sample was annealed 
at 118 K, the main SO 4 - band continued to increase, other weaker bands of SO 4 - 
(606 cm -1 and 980 cm -1 ) appeared, and the H 2 0 2 and S0 2 bands decreased signif- 
icantly. Fig. 2 shows the H 2 0 2 abundance in the sample as a function of time (50 K; 
t = 0) for annealing temperatures of 1 1 0, 1 1 8, and 1 30 K. In all cases the H 2 0 2 abun- 
dance began to decrease around 100 1 < and continued to fall as the sample was an- 
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Fig. 1. Evolution of a H 2 0 + S0 2 + H 2 0 2 sample (80:14:6) during heating to and 
annealing at 118 K. The bottom four spectra are labeled by the time elapsed after 
the sample had reached 1 1 8 K. 



Fig. 2. Evolution of a H 2 0 + S0 2 + H 2 0 2 sample (80:14:6) during heating to and 
isothermal annealing at 110, 118, and 130 K. All samples were deposited at 50 K 
and warmed at 1 K/min. Inset: Arrhenius plot for the H 2 0 + S0 2 + H 2 0 2 reaction. The 
line’s slope yields an activation energy of 53 kj/mole. 


nealed. For the samples held at 110 and 118 K, the H 2 0 2 abundance dropped by a 
factor of two and a factor of eight, respectively, after 1200 min, while only 
~20 min were needed at 130 1< for the H 2 0 2 abundance to drop below the noise le- 
vel. Blank experiments with H 2 0-H 2 0 2 showed that over these same time periods 
there was no H 2 0 2 loss, as expected from the negligible sublimation rate of H 2 0 2 at 
130 K. 

The inset of Fig. 2 shows an Arrhenius plot derived from six different annealing 
temperatures, which together yielded an activation energy of 53 ± 5 kj/mole. The 
rate constant, k, was derived from N = N 0 exp(-kt), since we observed that, to a first 
approximation, the H 2 0 2 abundance (N) dropped exponentially with time (t) at 
each temperature studied. In the higher temperature experiments, the reaction be- 
gan before we reached the annealing temperature, so in all cases we calculated k by 
equating it to ln( 2 )/ti/ 2 , where ti / 2 is the time at the annealing temperature needed 
for the H 2 0 2 abundance to drop to half of its original value. We note that as the 
abundance of H 2 0 2 has already dropped by a factor of two by the time we warm 
to ~125 K (see 130 K experiment in Fig. 2), we only used annealing temperatures 
below 122 I< to calculate the activation energy. The error that we give in the activa- 
tion energy is a combination of the difficulty of fitting the H 2 0 2 baseline, the narrow 
temperature range (12 K) for which we can obtain a value for ti/ 2 , and the fact that 
at the upper temperatures the reactions already have proceeded to some extent by 
the time the sample attained those temperatures. Future studies will focus on 
reducing the uncertainty in these measurements by extending the experiments to 
lower temperatures so that the activation energy can be more-confidently used be- 
low 1 00 K. 

4. Discussion 

4.1. Reaction chemistry 

Thermal reactions between S0 2 and H 2 0 2 have been studied extensively within 
the atmospheric-chemistry community due to the importance of such chemistry in 
removing S0 2 , in liquid, gas, and solid phases, from Earth’s atmosphere. The pri- 
mary reaction sequence is believed to begin with the formation of bisulfite 
(HSOJ) through (1) 

H 2 0 + S0 2 — > HSO 3 + H + (1) 

where the H + either rapidly or in concerted fashion attaches to a second H 2 0 mole- 
cule to form hydronium: 

H + + H 2 0 H 3 0 + (2) 

In our previous studies with H 2 0 + S0 2 ices, we found that these thermal reactions 
occur in the solid state even at temperatures below 100K (Loeffler and Hudson, 
2010). We now find, in addition, that in the presence of H 2 0 2 the bisulfite is oxidized 
to form the sulfate (SOj - ) ion (Martin and Damschen, 1981): 

H 2 O 2 + HSO 3 SO^ + H 3 0 + (3) 

Fig. 1 provides evidence that these same thermally-induced reactions occur in 
our H 2 0 + S0 2 + H 2 0 2 ices. Beginning with the 50-K spectrum and working toward 
those for higher temperatures, IR absorptions of both HSOJ (1035 cm -1 ) and SO 4 - 
(1092 cm -1 ) products can be seen starting at ~80 K during warming. Both bands 
continued to increase with temperature until the SO 4 - feature was one of the 
more-prominent remaining IR absorptions. The increase in abundances for the 
HSOJ and SOj - ions correlated with the decrease in absorption band areas of reac- 
tants S0 2 and H 2 0 2 (Fig. 1), supporting the simple reaction sequence already de- 
scribed. Comparing the H 2 0 2 loss for our different annealing temperatures, we 
see that not surprisingly the reaction proceeds the fastest at the highest tempera- 
ture. However, even at the lower temperatures, the H 2 0 2 abundance is still decreas- 
ing after the longest times studied (~1200 min), indicating that given enough time 
all of the H 2 0 2 will convert to SO 4 - via reaction (3). 

4.2. Icy satellite implications 

Infrared spectra of the surface of Europa, Ganymede, and Callisto all contain an 
absorption at 4.05 pm, which usually is assigned to S0 2 (Hansen and McCord, 2008; 
Hibbitts et al., 2000; McCord et al., 1998). In contrast, the 3.5-pm band of H 2 0 2 has 
been found only in Europa’s spectra (Carlson et al., 1999; Hansen and McCord, 
2008). Our experiments show that the abundance of H 2 0 2 in these surface ices de- 
pends on the presence of S0 2 , and that, conversely, sulfur dioxide’s abundance will 
be influenced by the presence of H 2 0 2 . Given these new laboratory results, we now 
turn to some of the previous observations of the jovian icy satellites. 

On Europa, although both H 2 0 2 and S0 2 have been detected they do not appear 
to be uniformly distributed. Hansen and McCord (2008) determined that the S0 2 
infrared feature was present on Europa’s trailing hemisphere. This agrees with pre- 
vious measurements of Europa’s ultraviolet reflectance spectrum (Domingue and 
Lane, 1998; Hendrix et al., 2011; Hendrix and Johnson, 2008; Lane et al., 1981), 
which showed a strong slope on the moon’s trailing side, attributed to S0 2 
(0.28 pm band center). It also is consistent with the sulfur implantation rate being 
an order of magnitude higher on Europa’s trailing hemisphere (Johnson et al., 2004). 
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Interestingly, the more-recent ultraviolet measurements indicate that S0 2 is not 
only more abundant on the trailing side but that it is actually absent on the leading 
side, where the majority of the H 2 0 2 detections have been made (Carlson et al., 
1999). Thus, based on our results it seems plausible that the distribution of S0 2 
and H 2 0 2 may be a result of excess S0 2 consuming any H 2 0 2 produced by radiolysis 
in the trailing hemisphere and excess H 2 0 2 consuming the S0 2 formed by implan- 
tation in the leading hemisphere. Finally, we point out that one recent study did de- 
tect both H 2 0 2 and S0 2 in a spectrum of an ice-rich region on Europa’s trailing side 
(Hansen and McCord, 2008). However, it is possible that these two molecules may 
be spatially separated, but sufficiently close to lie within the same pixel area of the 
detector. 

On Callisto and Ganymede the 3.5-pm absorption diagnostic of H 2 0 2 is absent, 
which could indicate that surficial S0 2 is widespread or that the H 2 0 2 abundance is 
simply much lower than on Europa (Hendrix et al., 1999; Hendrix and Johnson, 
2008). The former possibility is consistent with IR measurements of Ganymede 
and Callisto (McCord et al., 1998), which showed that each satellite’s spectrum con- 
tained the 4.05-pm absorption band attributed to S0 2 . However, we note that there 
has been some discussion as to whether another species, such as H 2 C0 3 or other 
carbonates (Johnson et al., 2004), could provide an adequate match for the 4.05- 
pm feature in Callisto’s spectrum, as this band does not seem to be correlated with 
magnetospheric bombardment (Hibbitts et al., 2000). Also, the most recent ultravi- 
olet measurements show little variation across Callisto’s surface (Hendrix and John- 
son, 2008). 

As the H 2 0 + S0 2 + H 2 0 2 reaction we studied occurs quickly at temperatures rel- 
evant to these icy satellites (ti /2 ~ 1 yr at 100 K and ~0.3 h at 120 K), we expect that 
it plays an important role in the evolution of the jovian icy satellites’ surface chem- 
istry and may explain observations of Europa related to the presence and distribu- 
tion of H 2 0 2 as well as the lack of H 2 0 2 on Ganymede and Callisto. If other 
molecules prove to be reactive with H 2 0 2 at these or at even lower temperatures, 
then similar thermal chemistry may explain why H 2 0 2 has not been detected on 
most Solar System icy bodies exposed to radiation. Future laboratory studies will 
focus on extending the reported measurements to other temperatures and concen- 
trations, identifying and quantifying other thermal reactions that may occur be- 
tween H 2 0 2 and other astrochemically-relevant molecules, and establishing 
whether the reaction products contain IR absorptions which could be detected by 
remote sensing. Finally, we note that the reactions we have described are not re- 
stricted to the surface of Europa and other worlds. To the extent that vertical trans- 
port of radiolytic products occurs (Greenburg, 2010), our reactions (l)-(3) also will 
take place beneath the IR-sensed layer of Europa, providing a thermally-driven 
source of SO 4 - . 
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